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D Abstract 

We study an extension of the MSSM by an anomalous abelian vector multiplet 

CO ■ 

and a Stiickelberg multiplet. The anomalies are cancelled by the Green-Schwarz 

^5 mechanism and the addition of Chern-Simons terms. The advantage of this choice 

over the standard one is that it allows for arbitrary values of the quantum numbers of 
the extra U(l). As a first step towards the study of hadron annihilations producing 
four leptons in the final state (a clean signal which might be studied at LHC) we 



then compute the decays Z' — ► Z07 and Z' — ► ZqZq. We find that the largest values 
of the decay rate is ~ 10~ 4 GeV, while the expected number of events per year at 
LHC is at most of the order of 10. 
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1 Introduction 

The Standard Model (SM) of particle physics has been confirmed to a great accuracy 
in many experiments. Despite the fact that the Higgs particle remains experimentally 
elusive, few scientists doubt that there will be major surprises in this direction. The 
whole scientific community, however, knows that the SM needs to be improved. First 
of all, neutrino oscillation experiments have exhibited the evidence for (tiny) neutrino 
masses, that have to be incorporated in (an extension of) the SM. Many ideas exist on 
how this can be achieved and more experimental precision tests will indicate which models 
are viable. Second, there are also several theoretical issues that make physicists believe 
that the SM is only an effective manifestation of a more Fundamental Theory. 

In approximately one year, the Large Hadron Collider (LHC) at CERN will start to 
operate at energies of order of 14 TeV in the center of mass. Apart from the search for 
the Higgs boson, it will probably give us some answers about the parameter space of the 
physics beyond the SM. Among the many issues that will be addressed, it is worth to 
mention: the search for supersymmetry, heavy quarks and the quark-gluon plasma, the 
existence of extra dimensions and the possible creation of tiny black holes. 

One of the most attractive scenario for physics beyond the SM is the existence of 
additional massive neutral gauge bosons [l]-[9]. They could be one of the first discov- 
eries at LHC if their mass is in the range of a few TeV. Many different models have 
been developed in the past in order to investigate this possibility. The mass could be 
acquired in a variety of ways: from Kaluza-Klein modes to a standard Higgs mechanism 
or even by adding an axionic field, 0, which couples to the abelian factors (Stiickelberg 
mechanism) [10, 11]. The latter is common to low energy effective field theories which 
appear anomalous. The anomaly cancellation is achieved by the Green-Schwarz mecha- 
nism with Stiickelberg terms accompanied by axion like couplings, <pFF, which ensure 
the consistency of these models [12, 13]. 



For example, in string theory anomalous £7(1) 's are very common. D-brane models 
contain several abelian factors, living on each stack of branes, and they are typically 
anomalous [14]- [29]. In the presence of these anomalous C/(l)'s, the Stiickelberg mixing 
with the axions cancels mixed anomalies 7 [16], and renders the "anomalous" gauge fields 
massive. The masses depend non-trivially on the internal volumes and on other moduli, 
allowing the physical masses of the anomalous U(l) gauge bosons to be much smaller 
than the string scale (even at a few TeV range) [12, 30]. However, it has been shown 
that axionic terms alone are not sufficient to cancel all anomalies. An important role is 
played by the so-called Generalized Chern-Simons terms (GCS) which are local gauge non- 
invariant terms. Indeed, these trilinear gauge bosons anomalous couplings are responsible 
for the cancellation of mixed anomalies between anomalous U(l)'s and non anomalous 
factors ensuring the consistency of the theory [31] -[35]. 

In this paper, we are interested in anomaly related Z' bosons in a non-renormalizable 
effective field theory. More precisely, we study an extension of the MSSM (see [36] for a 
review) by the addition of an abelian vector multiplet V^ ** and we assume that generically 
all MSSM particles are charged with respect to the new U(l). In order to gain in flexibility, 
our model is only string inspired: we do not commit to a specific brane model and this is 
why the charges are not fixed, even if the effective cut-off is related to the mass of the Z' . 
The extra vector multiplet generically is anomalous and consistency of the model requires 
an additional Stiickelberg multiplet S with the proper couplings as well as GCS terms. 
As a consequence, the anomalous abelian boson becomes massive and behaves like a Z 1 '. 
Moreover, in order to break supersymmetry, we add the usual soft breaking terms and 
the new terms coming from the fermionic sectors of V^ and S. 

Our model contains many new features: new D and F terms (which are coming from 
the axionic terms and not from the GCS, in accordance with [35], due to the fact that 
the GCS's contain only vector multiplets in antisymmetric form), new couplings and new 
mass contributions in comparison with the MSSM. Explicit formulae are provide for all 
these terms in component fields. 

Since the Higgs fields might be charged under the anomalous U(l), a combination of 
the Stiickelberg and the Higgs mechanism makes the anomalous U(l) massive. An axi- 
Goldstone combination is eaten by the neutral gauge bosons and no physical axi-Higgs 
is left contrary to other studies on anomaly related Z' [13] and similarly to the case of a 
non- anomalous related Z' [10, 11]. 

We explicitly show how the anomaly cancellation mechanism works in our model 
before and after breaking the gauge symmetry. Before gauge symmetry breaking, only 
SM fermions contribute to the triangle diagrams. After gauge symmetry breaking, all 



7 Irreducible anomalies are cancelled by the tadpole cancellation. 



fermions that become massive still contribute to the anomalous triangle diagrams. Their 
contribution is cancelled by new diagrams which involve the Nambu-Goldstone (NG) 
boson exchange. 

In order to explore some phenomenological implications of our setting, we then analyze 
the decays Z' — f Z 7 and Z' — ► ZqZq. We numerically compute the decay rates as 
functions of the arbitrary U(l) charges and the mass of the anomalous U(l) gauge boson. 
We find a non-trivial dependence on all these parameters, estimating that the region 
that gives the largest values is for Mz> ~ 4 TeV, where the decay rate Z' — ► Z 7 is of 
the order of 1CT 4 GeV. These decays are part of the processes in which two colliding 
protons lead to a four lepton final state [37]. The final state is very clean and possibly 
measurable at LHC. Assuming a degenerate mass spectrum for the sfermions of about 
500 GeV we also estimate Nz', the expected number of Z' produced per year. We find 
that Nz> falls off exponentially with Mz>, so we shall focus on the case Mz> ~ 1 TeV and 
the most favorite decay Z' — > ZqZq. We also estimate the number of decays for 1 year of 
integrated luminosity which turns out to be Nz'->z z ~ 10 in the most favourite region 
of parameters. In a future work we will push our program forward and study this signal 
with the aid of Monte Carlo methods [38] . 

The paper is organized as follows: in Section 2, we introduce the vector multiplet, V^°\ 
the Stiickelberg multiplet and we provide the axionic and GCS lagrangians in superfields 
and in components. We then discuss the anomaly cancellation both in the unbroken and 
in the broken phase. At the end of the Section, we add all possible soft-breaking terms. 
In Section 3, we describe the model set up. In particular, we discuss the kinetic mixing 
terms which are coming from the axionic lagrangian and the D and F terms, pointing out 
explicitly the new contributions. We comment on the superpotential and we compute the 
mass terms for all the particles, pointing out the differences from the canonical MSSM 
setup. Finally, in Section 4, we study some phenomenomogical implications of our model. 
We consider the case in which the Higgs fields are uncharged with respect to the U(l)' and 
compute the decay rates for the two processes Z' — ► Z 7 and Z' — ► Z Z which should be 
relevant for the computation of hadron annihilations into four leptons. In the appendices 
we report the technical details and discuss the general case in which also the Higgs fields 
transform under the anomalous U(l)'. 

2 Preliminaries 

In this section, we discuss how to extend the Minimal Supersymmetric Standard Model 
(MSSM) to accommodate an additional abelian vector multiplet V^ and how to cancel 
the anomalies with the Green-Schwarz mechanism. We assume that all the MSSM fields 



are charged under the additional vector multiplet V^\ with charges that are given in 
Table 1, where Qi, L; L are the left handed quarks and leptons respectively while U^, D\, E^ 
are the right handed up and down quarks and the electrically charged leptons. The 
superscript c stands for charge conjugation. The index % = 1,2,3 denotes the three 
different families. H u ^ are the two Higgs scalars. 
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Table 1: Charge assignment. 

Since our model is an extension of the MSSM, the gauge invariance of the superpo- 
tential, that contains the Yukawa couplings and a /i-term, put constraints on the above 
charges 



Qu c = 


'- ~Qq - Qh u 


Qd c = 


'- -Qq + Qh u 


Qe c = 


-Ql + Qh u 


Qu d = 


'- ~Qh u 



Thus, Qq, Ql and Qh u are free parameters of the model. 



2.1 Anomalies 



As it is well known, the MSSM is anomaly free. All the anomalies that involve only the 
SU(3), £77(2) and U(l)y factors vanish identically. However, triangles with £/(!)' in the 



external legs in general are potentially anomalous. These anomalies are 8 
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where / runs over the fermions in Table 1, Qf is the corresponding U(l)' charge, Yf 

p(o) 



is the hypercharge and TJ? , a = 2,3; fc a = 1, . . . , dirnG^ are the generators of the 



G^ = SU{2) and G^ = SU(3) algebras respectively. In our notation Tr [T-T^] = \8jk- 
All the remaining anomalies that involve U(l)'s vanish identically due to group theoretical 
arguments (see Chapter 22 of [39]). Using the charge constraints (1) we get 

-4 (0) = z{Qh u +3Qh u Q 2 l + C&-ZC&AQl + 6 Qq)} (7) 

-4 (1) = -\(ZQq + Ql) (8) 

-4 (2) = I^Qq + Ql) (9) 

A i3) = (10) 

-4 (4) = -GQhAWq + Ql) (11) 

Notice that the mixed anomaly between the anomalous U(l) and the SU(3) nonabelian 
factors .A® vanishes identically. 



2.1.1 Anomalous U(l)'s and the Stiickelberg mechanism 

Many models have been developed in the past where all the anomalies (7-11) vanish by 
constraining the charges Qf (see [1, 2] and references therein). On the contrary, in this 
paper we assume that the U(l)' is anomalous, i.e. (7)-(ll) do not vanish. Consistency 
of the model is achieved by the contribution of a Stiickelberg field S and its appropriate 



8 We are working in an effective field theory framework and we ignore troughout the paper all the 
gravitational effects. In particular, we do not consider the gravitational anomalies which, however, could 
be canceled by the Green-Schwarz mechanism. 



couplings to the anomalous U(l)'. The Stiickelberg lagrangian reads [40] 

A e 2 e 2 

J2 b i a) S Tr (W^WW) + bfS W {1) W w + h.c. 1 (12) 

.o=0 J 6(2 J 

where the index a = 0, . . . , 3 runs over the U(l)', U(1)y, SU{2) and SU(3) gauge groups 
respectively. The Stiickelberg multiplet is a chiral superfield 

S = s + iV26^ s + 2 F S - iOa^d^s + —6 2 6a^d^s - -6 2 6 2 n s (13) 

and transforms under the U(l)' as 

V (o) ^ 0°)+i(A-At) 

S -> S - Ai b 3 A (14) 

where 63 is a constant. The lowest component of S is a complex scalar field s = a + i<f). 
We assume that the real part a gets an expectation value by an effective potential of 
stringy or different origin and contributes to the coupling constants as 

\bt\a) (15) 



I6glr a lQgiT a 2 

where g a is the redefined coupling constant and the gauge factors r a take the values 
1, 1, 1/2, 1/2. The first line in (12) is gauge invariant and provides the kinetic terms and 
the axion-[/(l)' mixing. The second line is not gauge invariant and provides couplings 
that participate in the anomaly cancellation procedure. Notice that in (12) the sum over 
a omits the a = 3 case since there is no mixed anomaly between the U(l)' and the SU(3) 
factors as from eq.(10), i.e. b 2 — 0. The values of the other constants, b 2 ° , are fixed by 
the anomalies. 

At first sight our lagrangian (see Appendix B) may look not the most general pos- 
sible one. In particular, an explicit Fayet-Iliopoulos term ^V^ could be added. It is 
well known that in certain string-inspired models (see, e.g. [41]), an one-loop FI term is 
absent, even if Tr(Q) 7^ 0. This is in apparent conflict with the observation [42] that in 
field theory a quadratically divergent FI term is always generated at one loop. The solu- 
tion to this paradox is that in the low-energy lagrangian there should be a counterterm, 
which compensates precisely, i.e. both the divergent and the finite part of, the one loop 
contribution. We do not write explicitly this counterterm, since its exact expression is 
model and regularization dependent, but we implicitly assume that such a cancellation oc- 
curs. As mentioned before, also the terms responsible for the cancellation of gravitational 
anomalies are omitted. 



Expanding C ax i on in component fields, using the Wess-Zumino gauge and substituting 
a by its vev we get 
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+ /i.e. 



(16) 



where we omit terms which are coming from (a)W^W^ a \ since they are absorbed in the 
coupling constant redefinition (15). This mechanism cancels some mixed anomalies and 
in addition provides a mass term to the anomalous U(l). Therefore, the anomalous U(l) 
behaves almost like the usual Z' extensively studied in the past. 



2.1.2 Generalized Chern-Simons terms 

As it was pointed out in [33], the Stiickelberg mechanism is not sufficient to cancel all the 
anomalies. Mixed anomalies between anomalous and non-anomalous factors require an 
additional mechanism to ensure consistency of the model: non gauge invariant Generalized 
Chern-Simons terms (GCS) must be added. In our case, the GCS terms have the form 

[32] 



C GCS = -<U[{VWD°VW -V<®D a vW)wW + h.c.] ff2S 
+d 5 [(vWD a V<® - V<MD a vW) W$> + h.c] m 

+ l -V^D a V^D 2 {[D a V^\V^])+h.c. 



+ 
+ 



(17) 



The constants d 4 , d§ and d§ are fixed by the cancellation of the mixed anomalies. The 
GCS terms (17), expressed in component fields, are 



+d 6 e^ pa V^ 0) Tr 



Cgcs = -d 4 e^vPv,^F^ + d 5 e^V^V^F^ 

vS*FV-ivP[v?\vV] 

-d 4 (\W<t»\<®VW - \<®<T"\WVW + h.c.) 
+d 5 (\Wa»\WvW - \Wa"\WvW + h.c.) 
+d 6 Tr [A(VA( 2 V;> - \WaV\WvW + h.c] (18) 

These terms provide new trilinear couplings that distinguish these models from the Z' 
models studied in the past. 

2.2 Anomaly cancellation 

In the following, we illustrate the anomaly cancellation procedure both in the unbroken 
and broken phases by a specific example. We focus on the bosonic sector and the related 
diagrams, since their supersymmetric analogs are fixed by supersymmetry. The GS and 
GCS terms depend on unknown parameters which we fix by using the Ward identities. In 
theories with massive gauge bosons where the mass is acquired either by the Higgs or by 
the Stuckelberg mechanism, Ward identities have the following diagrammatic form [43] 



-ik" \ V»(k) >~JM\Pl +m v \ G v (k) H$1PI =0 (19) 




where V^ is the massive gauge field, Gy is the corresponding Higgs or Stuckelberg field 
(or a linear combination of them) and my is the coupling of the term V^d^Gy. The blob 
denotes all the 1PI diagrams. 

2.2.1 Anomaly cancellation in the symmetric phase 

In our model there are two extra states in the neutral fermionic sector, namely the axino 
and the primeino (see Section 3.7) which do not contribute to the fermionic loop. The 
remaining MSSM fermionic states are a bino, a wino and the two higgsinos. Both U(l)y 
and SU{2) gauginos do not contribute to the fermionic loop due to group theoretical 
arguments (see Section 28.1 of [44]). The higgsino eigenstates do not participate because 
the H u contribution is cancelled by the Hj one. This is due to the fact that each diagram 



(p + q) p Vi°\p + q) 



p" v {0) 




q u V<® 




Figure 1: The Ward identities for the amplitude V p (p + q) — ► V^ (p) V„ (q) in the 
unbroken phase include the GCS as well as the axionic couplings. The solid lines represent 
fermions and the wiggle lines are gauge fields. Dashed lines are scalars. Each depicted 
diagram also contains the exchange (//,p) <-> {v, q). 

is proportional to an odd product of charges and the two higgsinos have opposite charges 
(see Table 1 and the constraints (1)). Without loss of generality, we assume that the 
mixed anomaly between V^ and two V^ 1 -* is non vanishing, therefore from eq. (3) A^ = 
YlfQf(Yf) 2 7^ 0- I 11 order to cancel the anomaly, we have to satisfy the Ward identities 
which are shown, in diagrammatic form, in Fig. 1. The total fermionic triangle is given 
by 

A™ 



1 *-^ A (l> 

a^(p, q; o) = -- J2 QA Y ff t pAp, ?; o) = -^g- r^ip, q] o) 



(20) 



The superscript indices in the l.h.s. stand for the gauge groups of the vector fields involved 
in the process. T ppu (p, q; 0) can be parametrized as in (134). For a symmetric distribution 
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of the anomaly (see Appendix C.2), we have 



(p + qyA™(p,q\0) = ^§^WPV 



3 32tt 2 
1A™ 



Denoting by 



1 AW 
<f±™(p,q',0) = -^e^V (21) 



(OS) 1 ' = -2a£ ) e ltmt arq f> (22) 



the axion interaction vertex and by 

(GCS) p f l/ = 2d 5 e pufia (p-qy (23) 

the GCS coupling, the Ward identities in Fig. 1 correspond to 

(p + qy(A™(p,q;0) + (GCS)™)+2ib 3 (GS)l 1 l/ = 
p^A p f u (p, q -,0) + (GCS)f l ! u ) = 
<r(A o £,(p,q-,0) + (GCS)™) = (24) 

They fix the parameters 



., hj| __^« J^_ 

1287T 2 1927T 2 



bZ'ba = r d 5 = r (25) 

2 3 128tt 2 5 192tt 2 V ; 

In the same way, the cancellation of the remaining mixed anomalies gives 



A A(0) h w h A(2) h w h A(4) 



384vr 2 2 64tt 2 l ° 128tt 2 

d 4 = - J ^ d 6 = : — 7 (26) 

384vr 2 6 96tt 2 v ; 



It is worth noting that the GCS coefficients (£4,5,6 are fully determined in terms of the A's 
by the Ward identities, while the b 2 s dep 
related to the mass of the anomalous U(l). 



by the Ward identities, while the b 2 s depend only on the free parameter 63, which is 



2.2.2 Anomaly cancellation in the broken phase 

It is interesting to study the anomaly cancellation procedure in the broken phase. Fo- 
cusing again onto the non- vanishing A^ 7^ 0, the amplitudes that contribute to the 
cancellation of the anomaly are given in Fig. 2, where mo = Qh u \v\/2 and mi = \v\/A 
with \v\ = \/v^ + v\. In the broken phase, additional contributions coming from the NG 

11 
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Figure 2: The Ward identities for the amplitude V p (p + q) 
broken phase. 




boson exchange must be added. We denote by A mu (p,q;mf) the modified triangle di- 
agram where also massive fermions circulate in the loop and by (NG) ppu the triangle 
diagram with a NG boson on an external leg. Note that (GS) ppu and (GCS) ppu are the 
same as in the unbroken phase. The amplitude satisfies again the usual Ward identities 
(24). In order to clarify the mechanism, we will focus on a single Ward identity 



(p + ?)" (A^(P, q\ m f ) + (GCS)™) + 2ih(GS)H + zm (NG) 
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From now on the (p, q; mf) dependence will be explicit only when needed. Splitting A 
and (NG) terms into the sums over SM fermions and higgsinos we obtain 
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Table 2: Definition of £? n and tf GU , where iV c = 3 is the number of colours. 



Since we have 
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where the integral Jo is 
I (p,q;m f ) = - 



1 pl—x 

dx / dy 

i) Jo 



y(l — y)p 2 + x(l — x)q 2 + 2xyp ■ q — ml 



(31) 



and t°r 11 , t^ Gl1 are defined in Table 2, the Ward identity of the SM fermionic loop has a 
new contribution due to the masses of the fermions. Similarly, for the corresponding NG 
term we get 



im (NG) 



n I 

M*dsM 



48tt 



l - 2 y, [*r ' -j '«] wv 



(32) 



feSM 



Summing (30, 32), the massive contribution in the fermionic loop is exactly cancelled by 
the NG ones, giving 



[(p + ?)'A™ (p, g; m/) + imo(iVG)^] 



SM 



1 V^ + 011 c „«J 



96tt 
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(33) 



The contribution of the diagrams involving the higgsinos vanishes 

[(p + ?) p ASi(p,( Z ;m / )+»mo(JVG)a^ = ^ E Q/*/WpV = (34) 



96tt 2 



/eff„ 



Summing (33, 34) we get 



A™ 



[{p + g)"A^(p, g; m f ) + im^NG)*] = — e^pV = (p + qf A^(p, g; 0) (35) 

Thus the contribution to the Ward Identities of the triangle diagrams is exactly the same 
as in the unbroken phase. 
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2.3 Soft breaking terms 

The total soft breaking lagrangian can be written as 



with 



f r<MSSM , nnew ( r XCC\ 

J-Soft — '-'soft "T" '-'soft \ OD ) 



C/f M = -\ E (M a A«A« + ft.c.) - (mJ^Qj + KJJfUf + ml^D] 
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l Lij L i^j -^•"E i ,^; ± -j - '"h„\ 



+mjL t L\ + mVE c M + ml \h u \ 2 + m\\h d ? 



(a%QiU?K - aZQifyhd - a'fUE]^ + bh u h d + h.c^j (37) 



and 



C?t = ~\ (M \ (0) \ {0) + h.c.) - \ (M s ^ s + h.c.) (38) 



where A*- -* is the gaugino of the added U(l)' and tps is the axino. We allow a soft mass 
term for the axino since it couples only through GS interactions and not through Yukawa 
interactions [45]. Notice also that a mass term for the axion is not allowed since it 
transforms non trivially under the anomalous U(l)' gauge transformation (14). 

3 Model setup 

In this Section we analyze the effects of the additional terms on the rest of the lagrangian. 

3.1 Kinetic diagonalization of U(l)'s 

As we mentioned before, the Stiickelberg multiplet contains a complex scalar field whose 
real part gets an expectation value that modifies the coupling constant (15). Therefore, 
the second line in (12) contributes to the kinetic terms for the gauge fields and the term 
(a)W^W^ gives a kinetic mixing between the V^ and V^ gauge bosons. Redefining 
as usual V<® -► 2g V (0 \ V™ -► 2g{VW we get 



4 4 2 J 



(39) 

e 2 



with 5 = — Ab 2 gogi(cx). In order to diagonalize the kinetic terms, we use the matrix 

i: ) - ( 3, : ) ( x > 

14 



where Cs = 1/vl — S 2 and S$ = 5C$. Let us stress that in this case the mixing is a 
consequence of the anomaly cancellation procedure. Note that, since b 2 ~ b^ 1 ~ M~^ 0) 
(see eq. (26)), where M V (o) is the mass of the anomalous U(l) that we assume to be in 
the TeV range, this mixing is tiny and can be ignored for our purposes. 



3.2 D and F terms 

The additional fields give rise also to D and F terms. More precisely, D term contributions 
come from: (i) the kinetic terms of chiral multiplets and (ii) the axionic lagrangian, 
providing (12) 



)(0) 



a=0 



£>» = ^E^>£ } + E^i a >K T i: ) )^' +±9oh(a)D^ +8DW 

-0 a=0 

" 2 Z,(4) 

J2 9l 4 a) v/2^Tr (A«i?M) + 9o 9 X^s (X^D® + X^D^) + h.c. 

n — n v 



a=0 



(41) 



where a = 0, 1, 2, 3 denotes, as usual, the gauge group factors, Z{ are the lowest components 
of the i-th chiral multiplet (except the multiplet which contains the axion) and T k , 
k a — 1, . . . , dimG^, are the generators of the corresponding gauge groups, G^. Solving 
the equations of motion for the D's and substituting back we obtain 



£ 



D C 



Csgo J2 Qf\ z f\ 2 ~ S s9i ^2 Y f\ z f\ 



Cs4:g b 3 {a) 



+ 2v / 2~4%o fys {Cj\ c ) + h.c] + 2V2b^g 2 [fa (S 2 5 \ c - S 5 \ B ) + h.c] 
+ V2b^g gi [fa (C S X B - 2C 5 S 5 X C ) + h.c] } ( 12) 



£d t 



\\giY. Y f\ z f\ 2 + 2 ^ b 2 ] 9 2 i bPs (X B - S s Xc) + h.c] + 



+ V2b^go9i [ipsCsXc + h.c] 



£_D<2) 



EUW+W 



2faX^ ] + h 



,]} 



^EMctfV 



(43) 
(44) 
(45) 
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Similarly, the F term contributions are 

dW 



C 



E ( FfF i 

MSSM \ 



f&MSSM 
1 

2 



__ F / _ ^! F t 



\~F s Fl + ^ { F s 



1 
2 



^^Tr^A^H-^A^A^ 



a=0 



+ h.c. 



(46) 



where the first line is the standard MSSM F term contribution while the second line 
contains the new axionic terms. Solving the EOM, and rescaling V — ► IgV we get 



£ 



Fs 



^4 a) ^Tr(A (a) A (a) )+^rfA (1) A (0) 



(47) 



Eq. (47) can also be written in the basis (40), but we will not need this term in the 
following. 

We would like to mention that no D and F terms are coming from the GCS since they 
include only vector multiplets in an antisymmetric form. Our results are in accordance 
with [35]. 



3.3 Scalar potential 

As we have seen in the previous section, the additional F terms (47) do not give any 
contribution to the scalar potential. The D#, D^ 2 ^ and D® terms (see eq. (43), (44) and 
(45)) provide the usual contributions to the MSSM potential. The only new contribution 
comes from the first line of (42). Thus the scalar potential can be written as 



V 

Vd c 



Vmssm + Vd c 

1 



C s g J2Qf\ z f\ 2 - s s9iJ2 Y fM 
f f 



+ C s 4g b 3 (a) 



Solving the equations for the minima of the potential 

dV 



dz 







(48) 
(49) 



/ 



(50) 
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we get (zf) = for all the sfermions as in the MSSM case. Inserting back these vevs into 
(48) we get the following Higgs scalar potential 

V h = {\»\ 2 + m 2 hu + 4g 2 h{u)CsX 6 }(\h u \ 2 + \K\ 2 ) 

+ {H 2 + m 2 hd - Ag 2 Ma)C s X s ) (|^| 2 + \h~ d \ 2 ) 

+{i(^oX,) 2 + i(, 2 + , 2 )}(|^| 2 + |^| 2 -|^| 2 -|^| 2 ) 2 

+ {HKK - h°X) + h.c) + \g 2 2 \KK* + h° u h-/\ 2 (51) 

which can be brought to the same form of the MSSM potential, after the following redef- 
initions 

ml u + 4g 2 h(a)C s X s -»• m\ u 
m l d -^3lh(a)C 5 X 5 -> m\ A 

({9oXs? + \{9l + gl)}v 2 - m 2 z (52) 

where 

g X 5 = C s g Q Hu - -S s gi (53) 

At the minimum, we recover the MSSM result (/i+) = (/ij) = for the Higgs charged 
components. Defining {h®) = Vi/y/2 , t> 2 + v\ = v 2 and v u /vd = tan/3 we can still write 
the tree level conditions for the electroweak symmetry breaking as 

b 2 > (H 2 + 0(1/^ + (54) 

2b < 2|/i| 2 + m 2 hu + m 2 hd (55) 

in complete analogy with the MSSM case (using m's). 

3.4 Higgs sector 

It is worth noting that in our model there is no axi-higgs mixing. This is due to the fact 
that we do not consider scalar potential terms for the axion (on the contrary to [13]). 

After the electroweak symmetry breaking we have four gauge generators that are 
broken, so we have four longitudinal degrees of freedom. One of them is the axion, while 
the other three are the usual NG bosons coming from the Higgs sector. 

As it was mentioned above, the potential has the standard MSSM form, upon the 
redefinitions (52). The Higgs scalar fields consist of two complex S'C/(2)i-doublets, or 
eight real, scalar degrees of freedom. When the electroweak symmetry is broken, three 
of them are the would-be NG bosons G°, G ± . The remaining five Higgs scalar mass 
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eigenstates consist of two CP-even neutral scalars h° and H°, one CP-odd neutral scalar 
A and a charge +1 scalar H + as well as its charge conjugate H~ with charge — l. 9 The 
gauge-eigenstate fields can be expressed in terms of the mass eigenstate fields as 



(56) 

(57) 



where the orthogonal rotation matrices R a ,Rj3 ,R/3 ± are the same as in [36] Acting with 
these matrices on the gauge eigenstate fields we obtain the diagonal mass terms. Ex- 
panding around the minimum (56) one finds that (3$ = (3± = P, an d replacing the tilde 
parameters (52) we obtain the masses 





in 



A" 



m h o h q 



21 1 2 i 2 i 2 

W + m h u + m h d 



(5* 



=F 



m% 



1 



(g X 5 y + -(gf + g 2 2 ) 



+ I ( (g X 6 ) 2 + -(gl + gl) ) v 2 m\ sin 2 (2(3) 



in 



H± 



in 



A° 



+ m 



w 



in 



AO +92' 



(59) 
(60) 



and the mixing angles 



sin 2a 

sin 2(3 

tan 2a 

tan 2(3 



m\m + m 2 



"ho 



ii" 



2 2 



m 2 A0 +({g Q X s ) 2 + \{gl + gl))v 2 
m\ -({goX s ) 2 + \{gl + gl))v< 



(61) 



Notice that only the h° and H° masses get modified with respect to the MSSM, due to 
the additional anomalous U(l)'. 

3.5 Neutral Vectors 

There are two mass-sources for the gauge bosons: (i) the Stiickelberg mechanism and (ii) 
the Higgs mechanism. In this extension of the MSSM, the mass terms for the gauge fields 



3 We define G = G + * and H = H + * . Also, by convention, h° is lighter than H°. 
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are given by 



C 



"' ~~ 2 \ ** Bfl V ^ 



M z 






V. 



(2)A* 



(62) 



C^, 5^ are the lowest components of the vector multiplets Vc, Vb- The gauge boson mass 

matrix is 

gogi^-Xs - g g 2 ^X 5 \ 



(Ml 



M 2 



2 
2i> 2 

9iT 



~9i92- 



\ 



92 



2v A 



(63) 



16g 2 b 2 C 2 + g 2 (v 2 )X 2 and the lower dots denote the obvious terms under 



where M 2 < 

symmetrization. After diagonalization, we obtain the eigenstates 



A, 



j 0m 



Z'„ 



g 2 B fl + g 1 v£ ) 

VgJ 



9l 



9iB^ 



9\ 



9l 



+ goQ 



tiu 



M 



g\v 2 



1M 2 



C, + O[g 3 ,M, 



3 1 

y(o)J 



C 



9oQh u v 2 



2 /If 2 



9\B U 



n V {2) 



O[g 3 ,M } 



3 1 



and the corresponding masses 



Mi 



M 



z 



Ml 







- A (9l + 9l)v 2 -{Q Hu 



Kv»+9o 



(Q 



H u 



,2 (9J + gf) ffo^ 4 

4^ (0) 

4M 2 






647T 2 M, 



y(o)_ 



v' + Ok&M, 



y(o)J 



(64) 

(65) 
(66) 

(67) 
(68) 

(69) 



where M V (o) = 45 3 g is the mass parameter for the anomalous U(l) and it is assumed to 
be in the TeV range. Due to their complicated form, the eigenstates and eigenvalues of 
M 2 (63) are expressed as power expansions in g and l/M V (o) keeping only the leading 
terms. Higher terms are denoted by 0[g^, M~f 0) ]. 

The first eigenstate (64) corresponds to the photon and it is exact to all orders. It 
slightly differs from the usual MSSM expression due to the kinetic mixing between V^ 
and \/«. 

For the rest of the paper, we neglect the kinetic mixing contribution since they are 
higher loop effects which go beyond the scope of the present paper. Then the rotation 
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matrix from the hypercharge to the photon basis, up to 0[qq, M. 
( Z' \ 



3 1 



is 



z 



0/i 



/ vi 0) \ 



o l3 



\ A » J 



V, 



(1) 



v v$ 



(70) 



9oQh u v 2 



_ n y/A+al 



v(°) 
3i 



9()Qh u V 
~92-™3 

v(o) 

32 



v(o) 







32 



31 



/ v; (0) 



v; 



(i) 



Vsf+sf / 



V v£> 



where i, j = 0,1, 2. 



3.6 Sfermions 



In general, the contributions to the sfermion masses are coming from (i) the D and F terms 
in the superpotential and (ii) the soft-terms. However, in our case, the new contribution 
comes only from the Dq terms 



v d c 

' 111 '.1 o o 



C s goQn u + -S s gi 



1„ \ fvl-v 2 d 



4Csg h(a) \\ J2 ( C s9oQf ~ S s9l Y f ) \y f \ : 
f 



(71) 



where the y/ stand for all possible sfermions. 



3.7 Neutralinos 



■-^ym^ + h.c. 



(72) 



With respect to the MSSM, now we have two new fields: ips an d A^ - 1 . Thus, we have 

c l 

neutralino mass 

where 

(r) T = (1>s, Ac, A s , A< 2 \ h° d , hi) (73) 

The neutralino mass matrix M^ gets contributions from (i) the MSSM terms, (ii) the 

h — h — A^ -* couplings, (iii) the new soft-breaking terms £™f t , (iv) the Stiickelberg action 

and (v) the D terms. Finally, we obtain the symmetric matrix 

/ M s m sc m SB ^cftbf Av 2 



M 



A? 



. . . M CI + M.S] 


-MiSs 







Mi 




M 2 



-g v d X s 
_grvd 

2 
92Vd 
2 



\ 
9qv u X s 



givu 
2 

2 



/ 



(74) 
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(If 


"? 


«/° 


W f 


°f 


Ve, V(i, V T 





1/2 


1/2 


Ql 


Ql 


e, /i, r 


-1 


-1/2 + 2 sin 2 # w 


-1/2 


2Ql 





■u, c, £ 


2/3 


1/2 -4/3 sin 2 # w 


1/2 


2Qq 





d, s, b 


-1/3 


-1/2 + 2/3 sin 2 ^ 


-1/2 


2g Q 






Table 3: Couplings of the SM fermions with the neutral gauge bosons. 



where M 1; M 2 are the gaugino masses coming from the soft breaking terms (37), and 



m sc = V2I 2 (c 2 g 2 bf + S 2 glb? - C 5 S s9m bf) (g X s Av 2 + C s M v(0) (a)) 
+\ {-2S 5 glb^ + Csgogib?) 9l Av 2 + ^M, 



v(°) 



(75) 



m SB = V2{ (Cggogib^ - 25^6^) (g X s Av 2 + C 5 M y(0) (a)) + b? g\ At 



with Av 2 = v 2 — v 2 d . It is worth noting that the D terms and kinetic mixing terms are only 
higher order corrections and they can be neglected in the computations of the eigenvalues 
and eigenstates. 



4 Phenomenology 

In this Section we compute the amplitudes for the decays Z' — ► Z$ 7 and Z' — ► Z Z 10 
focusing for simplicity on the case Qh u — 0. In this case there is no mixing between 
the V^ and the other SM gauge fields therefore Z' = V^ (see (70)). Notice also that 
neutralino and chargino contributions to the fermionic triangles identically vanish, giving 
the same results, for what the decays of interest are concerned, of non-SUSY models. In 
Table 3 we list all the couplings of the SM fermions with the neutral gauge bosons where 
qj denote the electric charges, t>,° and a f ° are the vectorial and axial couplings with Z 
and vf and of are the vectorial and axial couplings with Z', respectively (see also (156)). 

4.1 Z' -> Z 7 

We compute all the relevant diagrams in the R^ gauge, thus removing the interaction 
vertex V^d^Gy that involves the massive gauge bosons and the Stiickelberg or NG boson. 
Therefore, the only diagrams that remain are the fermionic loop, the GCS vertex and a 



We would like to acknowledge T. Tomaras for discussions on this point. 
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Z'P(p + q) 



7"(?) 




+ Z' 




Not An. 



Figure 3: Diagrams for Z' — > Zq 7. 



not anomalous remnant contribution (Fig. 3). It is possible to show that the last blob- 
diagram, that involves several diagrams, is equal to zero. For the interested reader we 
give further details in Appendix C.3. 

The decay rate for the process is given by 

Pf 



r (z' - z 0l ) 



32tt 2 M|, 



A 



TOT 



l dn 



(76) 



where Atot is the total scalar amplitude and pf is the momentum of the outgoing vectors 
in the CM frame 



Pf 



Ml 



(77) 



The square of the total scalar amplitude is given by 

|^TOT| - Z^ e (A') e (A') Z^ e (A°) e (AO) 2^ e (A^) e (A^) A WMi"l ^2/^2 l' f 



A" 



where e are the polarizations of the gauge bosons, and A piiv is the Feynman amplitude 
of the process. The factor 1/3 comes from the average over the Z' helicity states. The 
polarizations obey to the following completeness relations 

h pl h P2 



e (A') e (A': 



E 

A" 

E 



£ (A°) e (AO) 



-7/ 



PlP'2 



-V\V2 



+ 



"(A'H-V) 

M\, 

ft, (A°) /i (A0) 



c Pl *M2 

fc (A^) fc (A^) 



-?7 



MiA*2 



(79) 

(80) 
(81) 



where (81) gives only the relevant part of the sum over helicities. Other terms are omitted 
since they give vanishing contributions to the decay. 
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The amplitude is given by the sum of the fermionic triangle A^f 07 plus the proper 
GCS vertex 

aZ'Z 01 _ AZ'Zo-y , t r pa\Z'Z^ 

^-ppv — ^ppv "T \ UUJ )pHV 

A p ^ = -\g,gz eY,vfaf q} V™{ PA -m } ) (82) 

/ 

where rY^, y (p, q; rrif) is given by (129). It is convenient to express the triangle amplitude 
by using the Rosenberg parametrization [46] 

A Zu 01 = ~ j^#o#z e(Aie[p, p, is, p] + A 2 t[q, p, is, p] + A 3 e[p, q, p, p]p u 

+A 4 e[p, q, p, p]q u + A 5 e[p, q, is, p]p^ + A 6 e[p, q, is, p)q^j (83) 

where 

A i = J2 vfafqfk for i = 3, . . . , 6 (84) 

/ 
J 3 , J 4 , J 5 and 16 are finite integrals (their explicit forms are given in (135)) and e[p, q, p, a] 
is defined after (134). A\ and A 2 are naively divergent by power counting and so they 
must be regularized. We compute them by using the Ward identities. In this way it is 
possible to express A\ and A 2 in terms of the finite integrals I 3 , I4, I5 and Iq. The GCS 
term has the following tensorial structure 

d z ' z ^(e\p,p,v,p\-e[q,p,is,p\) (85) 

so it can be absorbed by shifting the first two coefficients of the Rosenberg parametrization 
for the triangle. The resulting amplitude can be written as 

A fpu 01 = -^9ogz e(Ai€\p, p, is, p] + A 2 e[q, p, is, p] + A 3 e[p, q, p, p]p u 

+A A e[p, q, p, p]q u + A 5 e[p, q, is, p]p^ + A 6 e[p, q, is, p]g M J (86) 

The Ward identities (19) on the amplitude now read 

(p + qyA z p '^ + iM z ,(GS) z y = (87) 

P »A z p '^ + iM Zo (NG) z ^ = (88) 

^<f 07 = (89) 

where M z > = 4& 3 g and M Zo are the Z' and Z Q masses respectively. After some manipu- 
lations we obtain 

{p + q) p A z p '^ = l- 2 g Q g Zo e l -Y^v z 'afq f e[p,q,p,v} (90) 

/ 

P^ppu 0y = -^9o9z eJ2 v f' aZ f 0( lf m2 f Ioe[q,p,is,p] (91) 

/ 
f A *z<n = (92) 
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and inserting (86) into the above identities we get 

A 1 = (q 2 A 4 +p-qA 3 ) 

M = (p 2 A 5 + P -qA 6 + (NG) z '~< 

with 



(93) 



(94) 



where I is the integral given in (31). Substituting A ± , A 2 from (93) into the amplitude 
(86) and performing all the contractions we finally obtain 

I A z' Z(n] 2 ,2 M> - Ml) 2 (Ml + M%) 



96M| o M|,7r 4 



J2vf'a Z f q f ((h + h)M 2 Zo +m 2 f Io 



(95) 



4.2 Z' -> Z Z 

The computations are similar to the previous case so we point out only the differences 
with the other decay. Mutatis mutandis, the decay rate for the process is given in (76) 
with 

M z , 



Pf 



1- 



4M|o 
M 2 



(96) 



The square of the total scalar amplitude is given by 



v\ *u 2 \ "* Ml *P2 aZ'ZoZq a*Z'Z Z 



\ A 1 2 _ ^ \ "pi *P2 

|^TOT| _ 3 Z^ (A') (A') Z^ C (A°) C (A°) Z^ C (A°) C (A") ^Pipi^i ^P2M2^2 



e,t^e, 



e^e^l A 



a: 



(97) 



A' A" A 

where the amplitude A pilv is always the sum of the fermionic triangle and the (GCS) term. 
The contribution to the fermionic triangle is 

1 



A 



Z'ZoZ 
Piiv 



2 

;9o9z 



f 



EWr 



VAV , Z' Z Z rVW 
pp^ + ^/ ^/ «/ J- ppv 



+ V (a z 'v z °v z ° T AVV + a z V°a z ° r AA4 ^l 

~ / j \ u n v n u n L ppv ' u n u n u n L ppv J 



) 



(9* 



where n runs over all the neutrinos while the r pA „,'s are given by (128), (131), (129), (130). 
Using the fact that for the three neutrino families we have v z = a z and v z ° = a z ° we 
write the total amplitude (the sum of the triangles plus GCS terms) as 

1 



A 



Z'Z Zo 
Pfiu 



S7T" 



9o9z ( A i e \P, V, v , p] + A 2 e[q, p, is, p] + A 3 e[p, q, p, p]p p 



+A 4 e[p, q, p, p\q v + A 5 e[p, q, is, p]p^ + A 6 e\p, q, is, p]q^ 



(99) 
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with 

Ai = 2 J2 vfafvfh for z = 3, . . . , 6 (100) 

/ 

where v z = 2v z for neutrinos and v z = vf for the other fermions. The Ward identities 
now read 

( P + qyA z J u ° z ° + tM z ,(GS)Zf» = (101) 

p>*A* z ° z °+iM Zo (NG)% z ° = (102) 

q v A z ' z ° z °+iM Zo {NG)lf = (103) 

leading to 

(p + q) p A z m z ° z ° = -Lgvgl^fafvfe^q^^} (104) 

/ 

^ A z'z z = _A_ gog 2 o J2v z 'a Zo v Zo m 2 f I e[q,p,u,p] (105) 

/ 

q v A z ; z » z « = -^g g 2 Zo J2vfa Zo v z °m 2 f Ioe[q,p, P ,ri (106) 

/ 

From these equations we find the following values for A\ and A 2 

h = (q 2 A,+p-qA 3 -(NG) z ' z ^ (107) 

A 2 = (p 2 A 5 + p-qA 6 + (NG) z ' z ^ (108) 

with 

{NG fz Q = J- vfafvf m)h (109) 

/ 
where Jo is the integral given in (31). Substituting back into the amplitude and performing 
all the contractions we finally obtain 

2 



lA z'z z 0l 2 _ 2 4 (Ml - AMI) 



192M^ tH 



J2 Vf'a Zo vf° ( 2(/ 3 + h)M 2 Zo + m)h 

. / 



(110) 



4.3 Numerical Results 



In this Section we show some numerical computations for the two decay rates T(Z' — > Zq'j) 
and T(Z' — > Z Z ). They depend on the free parameters of the model, i.e. the charges 
Qq, Ql and the mass of the Z' . We assume that Qh u = and we choose go = 0.1. 
We show our results in Fig. 4-6 in the form of contour plots in the plane Qq,Ql for 
M Z i = 1,2 and 4 TeV. Our choices for g , Qq, Ql and M Z i are in agreement with the 
current experimental bounds [37]. 
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The darker shaded regions correspond to larger decay rates. The white region cor- 
responds to the value 10~ 6 GeV that can be considered as a rough lower limit for the 
detection of the corresponding process. It is worth noting that increasing Mz' the mean 
value of the decay rate of Z' — ► Z 7 grows while the one of Z' — ► Z Zo decreases. We 
would also like to mention that increasing Mz> the iso-decay rate contours in the plot 
rotate clockwise getting more and more parallel to the Qi-axis. This effect is due to the 
fact that the contribution of the triangle diagram with the top quark circulating inside the 
loop becomes the dominant contribution for high Mz>- In this case the decays strongly 
depend on the top quark charge Qq while the lepton charges Ql become irrelevant. Fi- 
nally, we find that the region that gives the largest values (of order of 10~ 4 GeV) of the 
decay Z' — ► Z 7 is for M z > ~ 4 TeV and for Qq ~ 3, Ql ~ —2. 

To estimate the number of the anomalous decays that can be observed at LHC we 
shall use the narrow width approximation, 

A^particics = N z > BR(Z' -> particles) (111) 

where Nz> = (Jz 1 Ct is the total number of Z', BK(Z' — > particles) is the branching ratio, 
C = 10 34 cm _2 s _1 the luminosity and t =1 year. Finally oz> is the Z' production cross 
section [5] 

^ = ^^ E [UMUM + UMMx 2 j\nz' -> ft ft), (112) 

where / gii?i are the quark qi (or antiquark q{) structure functions in the proton, and the 
momentum fractions are 

xi,2 = {M z >/V~s)e ±y . (113) 

To estimate a rough upper bound for the anomalous BR we assume that the sfermions 
will have an universal mass of about 500 GeV. We integrate numerically the PDFs using 
a Mathematica package [48]. In Fig. 7 we show the result for Nz> at \fs = 14 TeV. We 
can see that the number of the Z' produced falls off exponentially with Mz', so we shall 
focus on the case Mz> ~ 1 TeV and the most favorite decay Z' — ► ZqZq. In Fig. 8, we 
estimate the number of decays for 1 year of integrated luminosity which turns out to be 
N z >->z Zo ~ 10 f° r large values of the charges Ql and Qq. We will present a more detailed 
analysis in a forthcoming paper [38]. 
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Z' -> Z y (M z , = l TeV) 



Z' -> Z Z„ (M z . = l TeV) 




' 
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Figure 4: M z < = 1 TeV. 



Z' -> Zo T (Mz' = 2 TeV) 



Z' ^ Z Z a (M z , = 2 TeV) 




' 
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Figure 5: M z > = 2 TeV. 



(M z . = 4 TeF) 




Z' -> Z Z (M z . = 4 TeV 



" 



Qq 





■ 
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Figure 6: M z > = 4 TeV. 
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Figure 7: Number of Z' produced at LHC in 1 year for C = 10 34 cm 2 s x and y^i = 14 
TeV, in units of Q%, in function of the mass of the Z'. 



(Qhu=0, go=0.1, M z .=l TeV) 
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Figure 8: Number of Z' -»■ Z Z at LHC in 1 year for C = 10 34 cm _2 s _1 , y/s = 14 TeV 
and M z > = 1 TeV. 
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A Conventions 

We use the space-time metric rj^ = diag(+, — , — , — ) and the spinorial conventions 

C2i = e 12 = 1 ei2 = e 21 = -1 en = e 11 = e 22 = e 22 = (114) 

V> Q = e Q % ^ a = e a ^ $ & = e & % $« = e & $ $ (115) 

^x = i ,a x a ^X = ^aX a (116) 

The Dirac matrices are 

„ ( o^ \ , f ^ = (1, -a) , . 

7 M = where <^ ' 117 

^ Oj | a" = (1,<?) 

and we define 

T5= ("oi) (118) 

B Total lagrangian 

The lagrangian of the model contains several terms 

£ = £<2 + £>L + C-gauge + £# + ^W + C-axion + ^GCS + ^So/t (H9) 

where 

C Q = {Q\e v(3) e v( \^\Q Qi y^ Qi 

+ (U^e- vi3) e-l vm e Q -f V(0) U: + {Dtfe-^e^e*^ D {) (120) 



C L = (Lle v( \-^ { \^ v(0) L l + (E^e v(1) e Q -f V(0) E^ e2 _ 2 (121) 



£* = (Hte v(2) e* vW e^ v(0) H u + H y (2) e^ vW e^ v(0) H d ) (122) 

£„ = (-LTr (^(3)^(3)) + 1 Tr (W^WW) 

+ Y^ Tr (^ (1) ^ (1) ) + Y^-pTr (IT(%(^ + /,c. (123) 

C-w = {y t JQ l U^H u -y t jQ l D c J H d -y^L t E^H d + f iH u H d ) d2 + h.c. (124) 
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c„ 



e 2 e 2 



' 2 

J2 b ( 2 ] S Tr (^(%W) + bfS W® W^ 



0=0 



+ /i.c. 



0- 



(125) 



£, 



GC5 



-d 4 [(V (1) D a V^ - VWD°VW) W® + h.c.] e2S2 + 

+d 5 [(vwirvw - vWD a vW) w& + h.c.] 02§2 + 

■HhTv (V^D a V^ - V^D a vW) W™ + 



~V^D a V^D 2 ([D a V^,V^])+h.c. 



\e 2 e 2 



(126) 



£ 



Soft 



l - I J2 M a \ {a) \ {a) + h.c. I - l - (Msiisips + h.c.) 

m 2 Qij QiQ* + mltfituf + m 2 Di .D?Df 
+m 2 Li l l L) + ml..E?Ef + m 2 h Jh u \ 2 + m 2 hd \h d \ 2 
- {(4QiU C jK ~ a i !Q l D c j h d - ai J LiE]h d + bh u h d + h.c. 



(127) 



where Cq, Cl and Ch provide the kinetic terms and the gauge interactions of the mat- 
ter particles such as (s)quarks, (s)leptons, Higgs(ino)s; C ga uge contains the kinetic terms 
for the gauge supermultiplet; Cw is the usual MSSM superpotential; C axion provides the 
kinetic term of the Stiickelberg multiplet and its Green-Schwarz interactions used in the 
anomaly cancellation procedure; Cgcs contains the Generalized Chern Simons interac- 
tions giving trilinear gauge boson couplings needed to complete the anomaly cancellation 
procedure; finally, Csoft contains the usual soft breaking terms of the MSSM as well as 
the new terms for the primeino and the axino. 

Notice that in order to include the coupling constants in the gauge interactions we 
need to redefine them as shown in equation (15) and to substitute V — ► 2gV. 



C Amplitudes, Ward identities and Anomalies 

C.l Fermionic loop diagram 

In this Subsection we give some general properties of the fermionic triangle diagram of 
Fig. 9. Consider a case in which only a single fermion circulates in the loop and each 
coupling is either axial (A) or vectorial (V) with charge equal to minus one. The fermionic 
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(p + q) f 




Figure 9: The anomalous triangle diagram. 



triangles containing an odd number of axial couplings, denoted by AVV , VAV, VVA and 
AAA are 



r™ v (p,q;m f ) = J J^)i Tr 



1 1 1 

p f- ~ 4 ~~ m f v f-~ rn f M f- + i> ~ m f . 

+ (p <->■ q,n<-+ v) 

i0 ' 1 1 1 

Ipl 1 1"! 757^ IT—, I + 



ITLf f — 771 f If. + f — 771 f 



rZ A (P, q ;m f ) = /(gj?> 



+(p <-> q,p <->• v 



+(q<r+ -(p + q)^ U ^p) 

U ' \ 1 1 

P f-~i~ m f U f — rrif ** If + p 1 — rrif 

+ (p^ -(p + q),l*<->p) 

1 1 1 

757p^ 1 757*1 757m 



771 f f. — 771 f If + f — 771 f 



128) 



(129) 



(130) 



:i3i: 



These integrals are superficially divergent (by power counting) and thus there is an ambi- 
guity in their definition. The internal momentum 0, can, in fact, be arbitrarily shifted (see 
Section 6.2 of [47]) 

£ a -^£ a + ap a + (a-(3)q a (132) 



leading to 



r fL v (p, q, P\ m f) = r tL v (p, ?; m f) - ^ e P»™(p - qY 



7T Z 



(133) 



The amplitudes (128), (129), (130) and (131) can be written using the the Rosenberg 
parametrization [46] as 

r p ^(p,q;m f ) = 

— [i^p, q; m f ) e[p, \x, v, p] + I 2 (p, q; m f ) e[q, \x, v, p] + I 3 (p, q; m f ) e\p, q, \x, p\p v 

+ h(p, q; m/) e[p, q, 11, p\q v + I 5 (p, q; m f ) e[p, q, v, p\p^ + I 6 (p, q; m f ) e[p, q, v, p\q^ 

(134) 



31 



with e[p, q, p, a] = t PiVpa p il q v and where 

-1 pl-X 



xy 



pi pL—X 

h(p,q;m f ) = - dx dy— -^- — ., 

Jo Jo y{l-y)p 2 + x(l- x)q 2 + 2xyp-q-m z 



pi pi—x 

h(p,q;rrif) = dx dy 

Jo Jo 



f 

1 r l-x 



X[X — 1 



l y(l — y)p 2 + x(l — x)q 2 + 2xyp ■ q — ml 

h(p,q;m f ) = -h(q,p;m f ) 
I 6 (p,q;m f ) = -I 3 (p,q;m f ) (135) 

In terms of the Rosenberg parametrization the /3 dependence of (133) is contained only 
in I\ and Ii ( which are superficially divergent). However, using the Ward identities, 



(p + q) p TZ v (p,q^;m f ) = ± 



- + m 2 f Io(p,q;m f ) 



e\p,q,n,v\ 



P^^ u v {p,q,(3;m f ) = — ——e[q,p,v,p\ 



2 + (3 



<f v pnv (p,q,P;m f ) = — g-2-e[g,P,P,A*] ( 136 ) 

where Jo is defined in (31), it is possible to show that they can be expressed in terms of 

I 3 . . . I & as 

2 + 8 
lf vv (p,q,(3;m f ) = p ■ q I 3 (p, q) + q 2 7 4 (p, q) + —^~ 

I? vv (p,q,P;m f ) = -lf vv (q,p, fcm f ) (137) 

From now on we omit the explicit j3 dependence to get more compact formulae. 



C.2 Anomaly distribution and cancellation. 

In this Subsection we show that the sum of the triangle amplitude and of the GCS vertex 

are independent of /3. Since the anomaly is independent of the fermion masses we discuss 

only the unbroken phase, i.e. irif = 0. We consider the anomaly between V^ and two 

V^. The total fermionic triangle (the sum of AAA+AVV+VAV+VVA triangles) can be 

written as 

A^ 
A°il(p, q; 0) = -— r^fo q; 0) (138) 

where A^ is the anomaly (3) and r Piiv is defined in (134). The Ward identities for the 
fermionic triangle are 

(p + qYA° p ll = _/^ WP V 



A^ 
<fA^ = P + fl^eprtfpP ( 139 ) 
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For instance, /3 = —2/3 corresponds to a symmetric distribution of the anomaly. The 
gauge invariance of the theory is restored using (see Section (2.2.1)) 

(p + q y(A^l(p,q ] 0) + (GCS)f l ! u )+2th(GS)ll = 
p^A^(p,q;0) + (GCS) 011 ) = 
<f(±™(p,q;0) + {GCS)™) = (140) 

The last two identities imply 

A^ 2 + /? A^ 

(2 + 0) -^— - 2d 5 = =*► rf 5 = + P - (141) 

v HJ 128tt 2 5 5 2 128tt 2 v ; 

and the first identity becomes 

A^ 2 + 8 AM m m .A {1) 

It is then clear that different choices in the anomaly distribution affect only the GCS 
coefficient d 5 while the GS coefficient b 2 remains the same. This means that removing 
the Stiickelberg coupling by gauge fixing and computing the physical amplitude A + GCS, 
we get the same result and the same Ward identity. Consider the amplitude 

A° p j! u = A™ + (GCS)™ = A^l + 2d 5 e pum (p - q) a (143) 

The GCS terms can be reabsorbed by the following redefinitions 

0^> = O 7 ^'- 2 * < i44 » 
O 1 *') = (S^«) +a * < i45 > 

Imposing the p^ and (f identities (24) we get 

h(p,q) = p-qh(p,q) + q 2 h(p,q) 

I 2 (p,q) = -h(q, P ) (146) 

that relate I\ and I 2 to the other Zj's. We can define 

f w „ = — (Ae[p, fi, v, p] + ^e[g, p, i/, p] + 7 3 e[p, Q, P, p]p" 

+he[p, q, fi, p\(f + I 5 e\p, q, u, p]p» + I 6 (e[p, q, is, p]<f ) (147) 

so that the amplitude is 



All = Kll + {GCS)™ = -^-f 1 ^ (148) 
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v 
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VI 
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Figure 10: Non Anomalous diagrams for trilinear neutral gauge boson amplitudes. 

and obeys the following Ward identities 

A™ 



P'Alll = 
fA™ = 



e^pV = -2tb 3 (GS)H 



(149) 



This result does not depend on the scheme of the anomaly distribution. 

C.3 Treatment of non anomalous diagrams 

In this section we show that the non anomalous diagrams in Fig. 3 vanish. The diagrams 
we consider, reported in Fig. 10, have no specific assignment for the external legs, to keep 
the discussion as general as possible. All the factors which are not relevant for our aim are 
omitted and all the possible leg exchanges are understood. Finally, we use dimensional 
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regularization and the R^ gauge with £ = 1, in such a way that each diagram vanishes 
separately. 

A) The Scalar triangle loop is given by 

n A f v f dM (2l+p-q) p (2l-q) u (2l+p) p 



wvr " y 7 (27r) 2 ^[(/-g) 2 -m 2 ][/ 2 -m 2 ][(/ + p) 2 -m 2 l 

d ^i ( 2 l+p-q) p (2l-q) l/ (2l+p)„ 



(2tt) 2w [(/ - q) 2 - m 2 ] [I 2 - m 2 ] [(/ + p) 2 - m 2 ] 

7 (27r) 2 ^[(/-p) 2 -m 2 ][/ 2 -m 2 ][(/ + g) 2 -m 2 ] l ' 

Performing the change of variable l p — ► — /^ in the second integral, one gets 
vW„ ^ _ / <^l {2l+p-q) p {2l-q) u {2l+p) p 



D ^p^ q > J (27r) 2w [(/-(/) 2 -//r][/ 2 -/;r][(/ + //) 2 -//; 2 l 

j2w 



/• ^ (2f + p- g ) p (2f+p) /i (2Z- g ) l/ _ 

7 (27r) 2 ^[(/-g) 2 -m 2 ][/ 2 -m 2 ][(/+j9) 2 -m 2 ] l J 



B) The "Scalar bubble loop" is given by 

d 2uJ l (2l+p + q) p r] pu 



D B pup {p,q) = -2 
= -2 



(2tt) 2w [I 2 - m 2 } [(l+p + q) 2 - m 2 } 

d 2uJ l (l+P + q) P V»v 

(2tt) 2w [I 2 - m 2 ] [(l+p + q) 2 - m 2 } 



J (2ir) 2 "[l 2 -m 2 }[(l+p + q) 2 -m 2 ] K ' 

Performing the change of variable I — ► — I — p — q in the second integral one gets 

d 2uJ l {l + p + q) p Vnv 



D*„ p (p,q) = -2 J 



(2ir) 2 " [I 2 - m 2 ] [(l + p + q) 2 - m 2 } 



+2 [ *?* ( l +P + l)p^ = o fi53) 

7 (27r) 2 -[(/+p + g) 2 -m 2 ][/ 2 -m 2 ] y J 

C) Since the ghost interact with neutral vectors only through the third component of 
SU(2), the Ghost triangle loop is proportional to 

^3bc^3cd^3db = — $bd e 3db = (154) 

The other diagrams in Fig. 10 can also be shown to vanish after manipulations similar to 
the ones used in (151), (153), (154). 



D Decay rates. General case 

In this Section we compute the amplitudes for the decays Z' — *■ Z 7 and Z' — ► Z Q Z in 
the general case Qh u 7^ 0, still neglecting the effects coming from the kinetic mixing. We 
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work in the limit 

g a v u , d « /i, M , M u M 2 , M s , M v(0) (155) 

in which msc ~ M V (o), msB ~ (see (75), (40), (53)). Hence, (74) takes the same form 
as in the symmetric phase in which neutralinos and charginos do not contribute to the 
anomaly (see Section 2.2.1). In the limit (155) an extension of the standard model by an 
extra U(l) and our SUSY model give the same results for what the decays of interest are 
concerned. 

We define the Dirac fermions \I// = I ) where Jl{r) are all the left (right) Weyl 

\f R J 
fermions in the model. The SM fermion interaction terms with the neutral gauge bosons 



arc 



C% = J&% = ~gz>*f'f(vj'-afT 5 )*fZl 

C-t = Jk Z ^ = -l9z o yfY{vf -a z f ° l5 )y f Z , 

£l nt = XA^-eqfVffVfAr (156) 



where 



< = <& + <& <$ = <%-<& 

vf = Qf L + Qf R af = Qf L - Ql 

Qf = Qf L = Qfn (157) 

The left and right charges are defined in the following way 

gz>Qf' L = g2T 3 2 + giY fL 01 +goQf L (158) 

gz'Qf R = giY fR O 01 +g Q fR (159) 

9z Qf° L = g2T 3 O 12 + g 1 Y fL O 11 +g Qf L O 10 (160) 

gz Qf R = 9iYf R O u +goQf R 10 (161) 

eQ fL = g2n0 22 + gi Y fL 21 = 9l Y fR 21 = eQ fR (162) 

(2) 

where Oij is given in (70) and T% is the eigenvalue of T 3 . 



D.l Z' -> Z 7 

The amplitude is given by the sum of the fermionic triangle Af t J : ° 7 plus the proper GCS 

vertex 

aZ'Z 01 _ \Z'Z ~f , (rr q\Z'Z -i 

*ZS" = -\9z'gz eJ2(vfaf qf r^ + afv^ qf rf^) (163) 

/ 
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The resulting amplitude can be written as 
aZ'z 01 = ' 



gz'9z e[Aie\p, p, v, p] + A 2 e[q, p, v, p] + A 3 e\p, q, //, p]p" 



with 



4tt 2 - 
+A 4 e[p, q, p, p\q v + A 5 e[p, q, v, p\f + A 6 e[p, q, is, p]q> 



Ai = J2 ( v f' a ? + a f v f°) If 1 * for i = 3, . . . , 6 



(164) 



(165) 



and the integrals ij given in (135). Ai and A 2 are the new coefficients with the GCS 
absorbed similarly to (146). 

The Ward identities (19) for the amplitude now read 

(p + q)>A%*" + iM z , [(GS)%? + (NG) %>] = (166) 

p»A z p '^ + iM Zo [(GS) z ^ + (NG) z ^] = (167) 

q v A z '^ = (168) 

where M z > and Mz are the Z' and Z masses respectively. In both (166) and (167) we 
have a (GS) and a {NG) contribution due to the two Goldstone bosons which are a linear 
combination of the axion and G°. We use (167) and (168) to fix A\ and A 2 while (166) is 
automatically satisfied. Contracting with p^ we get 



Hy^Z'Zo-y 



pfiu 



4g gl R**" b 2 % 3 + 2g gj B**" bfh + 2g 2 o9l R**" bf% 



+ -j-^ 9z ' 9z ° e ^2 v f' °/° q f m / J ° r e ^'^' z/ '/ ) ] 



4tt 2 

where Jo is the integral given in (31). The solution for A\ and A 2 is 

Ai = (g 2 A 4 + p • gA 3 ) 

i 2 = (p 2 A 5 + p ■ qA 6 ) + (GS) z '^ + (NG) z 'i 

with 

(NGf- 1 



(GS) 



Z'l 



Y,vfafq f m)l 

f 

327T 2 



+ 

(169) 

(170) 
(171) 

(172) 



9z'9z e 
The rotation factors are 



2 r>Z'Z Q1 1,(1) 



2 d^'^07 1,(2) I 



4^ ^oi 07 & 2 % + 2» 2 ^202 ° 7 &T'&3 + 2^i ^ooi 7 65 ^ 



,Z'Z 07 l(4). 



(173) 



7? 
R, 



Z'Z 01 
101 

Z'Z al 
202 

Z'Z 01 
001 



Oq\Oiq0 2 i 
Oq 2 0\q0 22 

o 10 o 21 



(174) 
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with Oij given by (70). Substituting A\, A 2 into the amplitude (86) and performing all 
the contractions we finally obtain 



I^tot|^z 07 - 9z>9z e QQMlMlir* 



X 



J2 If (vfaf + afvf) (J 3 + h)Ml + (GS) Z '^ + (NG) Z ^ 
L / 



(175) 



D.2 Z' ->• Z n Z 



^0 



The contribution to the fermionic triangle is 



/\Z'Z Z 



-g9z>9 Zo 



E 

/ 



7; Z V°7; Z ° T yAy + v z 'v Zo a Zo T VVA 4- fl Z V°7; Zo T AVV + 
v f a f v f l p^ + v f v f a .f l pi*> + a f v f v f *- pixu + 



+ «f«?«frr) 



(176) 



where the r p/UI/ 's are given by (128), (131), (129), (130). We write the total amplitude 
(the sum of the triangles plus GCS terms) as 



A 



Z'ZoZo 
ppv 



"^9z'9 Zo 



Ae\p, p, v, p] + A 2 e[q, p, v, p] + A 3 e[p, q, p, p\p h 



with 



where 



+A 4 e[p, q, p, p\q v + A 5 e[p, q, v, p}p^ + A 6 e[p, q, v, p\q^ 



A i = J2 t f Z ° Z ° I i fori = 3,..., 6 

/ 



t 



Z'Z Zo 



afvfvf + 2vfafvf + afafaf) 
and the integrals Ii are given in (135). The Ward identities now read 

(p + q) p A z ; z » z « + iM z> [(GS) z f° + (NG) z f°] 



(177) 
(178) 

(179) 



^A z ' ZoZo + iM 2 



P " 'ppv 



(GS) 



Z'Z 
pv 



(NG) 



Z'Z 
pv 



= 


(180) 


= 


(181) 


= 


(182) 



q v A z ' z ^ + iM Zo \(GS) z f + (NG) z f 



where Mz> and Mz are the Z' and Zq masses respectively. In (180)- (182) the (GS) and 
(NG) terms are present for the same reason as in the preceding Subsection. We use (181) 
and (182) to fix A\ and A 2 while (180) is automatically satisfied. 
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Contracting with p^ and q v we get 



v^ A z ' z ° z ° 



piiv 



\g\ R z m z ° z ° 4 0) 6 3 + 4 M 2 Rfoi° Z ° b { 2 ] h+ 



+ 



+2g gi R^° b^% + 2g£ gi i&f °'° 6^6 3 j 



183) 



1/ aZ'Z Z 



,3 d^'^o^o 1,(0) 



2 d^'^o^o l(1)i 



4^ i? V°"° ^6s + 4flb^ i?f f°"° #%+ 



H-2^ 2 i? 2 X ? oZo ft + 2<? V i?Sf oZo 6^63 



+ 



^9z'g 2 zoY,( v f' a ? v f° + §«f a ?°/°) m ? /o j e t^'^] 



;i84) 



where Jo is the integral given in (31). The solution for A\ and A 2 is 

A 1 = (q 2 A 4 + p ■ qA 3 ) - [(GS) Z ' Z ° + (NG) Z ' Z ° 
A 2 = (p 2 A 5 + p ■ qA 6 ) + (GS) z ' Zo + (NG) z ' z ° 



with 



(NG) Z ' Z ° 
(GS) Z ' Z ° 



E ( v ? a >f° + \<*>f°) m Po 



(185) 
(186) 

(187) 



/ 

64tt 2 
9z'9 2 Zo 



4tf} R Z o ZoZo bfh + ±9o9i Rfo ZoZo b { 2 % 3 + 



■29o9 2 2 Rlo2° Zo b?\ + 2g 2 9l Rg z ° z ° bf\ 



The rotation factors are 



RZ'ZoZo = 00 



L ooo 



io^io 



R z'z z = 



L 101 



01^10^11 



Rono ' — L^02^'10^'12 



L 202 



i2gfoZo = o w O 11 + O 01 O w O 



10 



;i89) 



with Ojj given by (70). Substituting back into the amplitude and performing all the 
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contractions we finally obtain 



12 2 4 ( M 1'- 4M |,) 

I^tot|z'z Zo = 9z'9zo — ig2M 2 4 — x ( 19 °) 

2 

J^*f Z ° Z °(h + h)M 2 Zo + (GS) Z ' Z ° + (NG) Z ' Z ° 
f 

(igr 
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